Moulding flexural waves in elastic plates lying atop a Faqir's bed of nails 
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Platonic crystals (PCs) are the elastic plate analogue of the photonic crystals widely used in 
optics, and are thin structured elastic plates along which flexural waves cannot propagate within 
certain stop band frequency intervals. The practical importance of PCs is twofold: these can be 
used either in the design of microstructured acoustic metamaterials or as an approximate model 
for surface elastic waves propagating in meter scale seismic metamaterials. Here, we make use of 
the band spectrum of PCs created with very small clamped holes, the nails of the title, to achieve 
surface wave reflectors at very large wavelengths, a flat lens, an endoscope, a directive antenna near 
stop band frequencies and cloaking from Dirac cones. The point pinned, Faqir, plate is particularly 
appealing as there is an exact dispersion relation available so the origin of these phenomena can be 
explained and interpreted using Fourier series and high frequency homogenization. 
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There has been much interest over the past 20 years 
in the analysis of elastic waves in thin plates in the con- 
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tinuum mechanics community 
an interest renewed in the metamaterial community with 
the theoretical proposals Hllli 
experimental validation [fj, [32 
and negatively refracting flat lenses for flexural waves. 

One of the attractions of platonics is that much of the 
physics of photonic crystals can be translated into pla- 
tonic crystals (PCs). There are mathematical subtleties 
in the analysis, and numerics, of the scattering of flexural 
waves [29I ] owing to the fourth-order derivatives, versus 
the usual second-order derivatives for the wave equation 
of optics, involved in the governing equations; even the 
waves within a perfect plate have differences from those 
of the wave equation as they are not dispersionless. 

There is a long history of wave propagation along peri- 
odically supported plates motivated by ribbed structures 
in underwater acoustics, (25|; the periodically pinned 
plate [24} being considered in detail. Movchan et al. 
[26J also give a complementary treatment for the flexural 
plate with doubly periodic domains of circular holes us- 
ing lattice sums and multipoles with a dilute limit giving 
approximate dispersion relations in closed form. Fourier 
series expansions introduced in 

[Hill and applied in @ 
also allow for a highly accurate analytic expression for 
the exact dispersion equation for pinned plates. We fol- 
low this latter route, which combined with the recently 
developed high frequency homogenization theory (HFH) 
uncovers the physics of pinned Faqir 34{ plates; a 
striking illustration (lensing via all-angle-negative refrac- 
tion) is shown in Fig. [IJa) , which contrasts with a shield- 
ing effect at a lower frequency in Fig. [T^b). The latter 
obviously has potential application in seismic metamate- 
rials [8[ for anti-earthquake systems. 

Much has been said about control of light [3l| , sound, 
water, or shear (SH) waves [9( using the rich behaviour 
encapsulated by the dispersion curves of bi-periodic 
structures, modelled by a Helmholtz equation, up to mi- 




Figure 1: (a) A forcing generates a bending wave totally re- 
flected by the array (with 22 clamped holes), throughout the 
zero-frequency stop band of Fig. [2] (here Q = 2). From the 
four-fold symmetry of the array and linearity of Eq[T] the de- 
fect at the center of the array (one missing hole) is shielded 
from bending waves generated by a forcing, (b) A bend- 
ing wave excited by a point forcing of normalized frequency 
Q = 2.6 is focused through an array (pitch 2) of 27 clamped 
holes (radius 0.01) tilted through an angle 7r/4. (c) Focus- 
ing through an array of 304 clamped holes of radius 0.01 in a 
regular square orientation at the same frequency of = 2.6. 
(d) Represents the focusing effects produced by an effective 
material governed by equation ((6| with T11T22 < yielding 
characteristic type of solutions. 



nor changes in the normalization of material parameters, 
and choice of boundary conditions (e.g. Dirichlet or Neu- 
mann for clamped or freely vibrating inclusions in the 
context of SH waves). However, when one moves into 
the area of elastic waves, governed by Navier equations, 
it is no longer possible to reduce the analysis to a sin- 
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gle scalar partial differential equation (PDE), as shear 
and pressure waves do couple at boundaries. There is 
nevertheless, the simplified framework of the Kirchhoff- 
Love plate theory [16j that allows for bending moments 
and transverse shear forces to be taken into account via 
a fourth-order PDE for the out-of-plane plate displace- 
ment field. This plate theory is a natural extension of the 
Helmholtz equation to a generic model for flexural wave 
propagation through any spatially varying thin elastic 
medium. It offers a very convenient mathematical model 
for any physicist wishing to grasp (some of) the physics 
of PCs using earlier knowledge in photonic or phononic 
crystals. However, while the Helmholtz equation can, 
with appropriate notational and linguistic changes, hold 
for acoustic, electromagnetic, water or out-of-plane elas- 
tic waves and so encompasses many possible applications, 
the Kirchhoff-Lovc plate theory is dedicated to the anal- 
ysis of flexural waves. For instance, it says little about 
propagation of in-plane elastic waves in platonic crystals. 
Bearing this restriction in mind, we solve 
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16l 22 1 for u(xi,x 2 ) on the square cell —1 < Xi,x 2 < 1. 
Here, n 2 = 12(1 - v 2 ) pu? / '{Eh 2 ) , where p, h, E, v are 
density, thickness, Young's modulus and Poisson's ra- 
tio of the plate, respectively, and uj is the angular wave 
frequency; the plate contains an array of very small in- 
clusions, with Dirichlet (clamped) boundary conditions, 
which is a good model for rigid nails of the Faqir plate. 

In the context of optics, ([1]) reduces to a second order 
PDE where only the first factor of ([1]) appears, in which 
case Dirichlet boundary conditions are a good model for 
infinite conducting thin wires at GHz frequencies, with 
the unknown u of such a Helmholtz equation being the 
longitudinal component of the electric field E z in TM po- 
larization and the spectral parameter f2 is associated with 
u) 2 e{x\ , X2 ) /c 2 whereby ui is the electromagnetic wave fre- 
quency, e is the relative permittivity and c is the speed 
of light in a vacuum. Such a model has been known since 
the mid 90s to lead to a zero frequency s top band asso- 



ciated with very low frequency plasmons [281 . [3 



For waves propagating through an infinite, perfect, 
doubly periodic medium, one can invoke Bloch's theo- 
rem [3, |2l| and simply consider a square cell with quasi- 
periodic conditions applied to the edges: 



u(l,x 2 ) = e 2lKl u(-l,x 2 ), u Xl (l,x 2 ) = e 2lKl u Xl (-l,x 2 ), 

(2) 



u(xi, 1) 



= e 2lK2 u(xi,-l), u X2 (xi,l) 



= e 2iK ^ 



2 u X2 (xi, -1), 

. . (3) 

with the Bloch wave-vector k — (ki,K2) characterizing 
the phase-shift going from one cell to the next. Also 
u = at the origin with continuity of u,u Xl ,u X2 there 



also. The exact solution is readily found 0, as 

^ exp(-«7rN • x) 

w(x) = exp^K-x) y 



[(Ki-Trni^ + Osa-Trna) 2 ] 2 -^ 2 ' 

(4) 

and enforcing the condition at the origin gives the dis- 
persion relation D(k, il) = where 



D( K ,n)= J2 



i 



[(71711 - Hi) 2 + (nn 2 - K 2 ) 2 



n 2 



As is well-known in solid state physics [7| only a lim- 
ited range of wavenumbers need normally be considered, 
namely the wavenumbers along the right-angled trian- 
gle TXM shown in the irreducible Brillouin zone inset 
in Fig. [5] There are however some exceptions, where 
it is possible to miss important details such as the stop 
band minima/maxima not occuring at the edges of the 
Brillouin zone @, 0, [li| and flat bands leading to strong 
anisotropy and slow light [l(| . 
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Figure 2: Bloch dispersion curves around the edges of the 
Brillouin zone TXM (inset) for an array of clamped holes 
(radius 0.01) of pitch 2. One notes the stop band for Q. 6 
[0,7r 2 /4] and three Dirac cones at frequencies Q. — ty 2 /2 (X 
point), ty 2 and 2tt 2 (F point). Solid curves are computed with 
FEM, and dashed curves are from HFH. 



The dispersion diagram is shown in Fig. [2] the singu- 
larities of the summand in Eq. ([5|) correspond to solu- 
tions within a square cell (without the nail) satisfying the 
Bloch conditions at the edges, in some cases these singu- 
lar solutions also satisfy the conditions at the support and 
are therefore true solutions to the problem. Solid lines 
in Fig. [2] are the dispersion curves double checked by 
the finite element method (FEM) using the commercial 
package Comsol, notable features are the zero-frequency 
stop-band and also crossings of branches at the edges 
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Figure 3: A bending wave excited by a forcing of normalized 
frequency 0, = 3.7952 (that of vanishing group velocity at 
r point in Fig. [2]) inside an array of pinned holes (pitch 2, 
radius 0.01) tilted through an angle tt/4 gives rise to respec- 
tively four (a) and two (b) highly-directed beams outside 24 
holes making a square (a) and 27 holes making a rectangle 
(b). When tilted through an angle 7r/4 (c) the same array 
with 27 holes acts has an endoscope for a source at frequency 
Q — 2.4674. The endoscope (d) and antenna (e) effects are 
simulations of the continuum PDEs ((6| generated by HFH 
with respective (Xii,T 2 2) coefficients of (-9.649,4.7149) at 
point M and (6.2524, 6.2524) at point F, where the effective 
media are highlighted by white lines. 

of the Brillouin zone. Importantly, HFH [j| captures 
the physics of this band diagram, with the asymptotic 
theory giving dashed curves that fit very well with the 
FEM (solid curves) near the band-gap edges which are 
the zero-group velocity and Dirac points. The multiple 
scales used for HFH are a short-scale = Xijl and a 
long-scale Xi — Xi / L for i = 1 , 2 where I and L represent 
respectively the characteristic small scale (half length of 
a cell) and the long scale. A small parameter is formed 
as e = l/L and an expansion of SI 2 = S1q + e£l\ + . . . 
and u = uq + tu\ + . . . is posed with respect to this pa- 
rameter. It turns out that the leading order term for u 
is wo(£,X) = /o(X)£/o(£) where the function /o repre- 
senting an envelope of the solution u is obtained by an 
equation only on the long scale as, 

/..../'..A A - ^/o = 0. (6) 

where Ty 's are integrated quantities of the leading order 
and first order short scale solutions with = for i =/= j 
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Figure 4: A plane bending wave of normalized frequency fl = 
9.7 (just below the first Dirac cone at T point in Fig. [5]) 
incident from the top on an array of pinned holes (pitch 2, 
radius 0.01) undergoes considerably less scattering (a) than by 
a clamped obstacle on its own (b). Again the asmptotic PDEs 
in equation (£9|) capture the essence of physics of cloaking (c) 
and scattering (d) by the same clamped obstacle. 

in the present illustrations and f^2 is simply the second 
order term of the Q expansion. As in |5| , assuming Bloch 
waves, the asymptotic dispersion relation for reads, 

Ti ■ 

0, ~ fio - (7) 

where Hi = Ki — di and di = 0, — 7r/2,7r/2 depending 
on the location in the Brillouin zone about which the 
asymptotic expansion originates, and fio is the standing 
wave frequency at the Brillouin zone edge ; these asymp- 
totics give the dashed curves in Fig. [5J For the case of 
multiple modes originating from the same point, as for 
example the Dirac points, Eq. ([7]) is no longer valid and 

(i) 

one obtains three coupled equations for /q that simplify 
to three similar anisotropic equations like ((6|). Note that 
Eq. ([6]) is of second order and not fourth. The fourth 
order plate equation yields two types of solutions namely 
propagating and exponentially decaying. Indeed in the 
long scale, solutions follow a second order PDE and the 
remaining information of the fourth order problem is en- 
closed in the Tjy integrated quantities. 

We give three potential applications of the platonic 
crystal described by fig. on top of the slab lens shown 
in fig. [T]as well as comparison solutions constructed from 
the continuum long-scale HFH theory; the coefficients Ty 
capture the potential strong anisotropy of the medium. 

Using the vanishing group velocity of the second dis- 
persion curve in fig. [5]in the neighborhood of the T point, 
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one can design a strongly directive antenna, as shown in 
fig. [3][a,b) and (e). The effective medium for the directive 
antennas is governed by Eq. © with T xl = T 2 2 = 6.25 24. 
In Fig. (3Je) the group velocity is calculated by differen- 
tiating Eq. ([7} with respect to Ki to obtain, 



0., 



Ill 

~n 



Ki, 



(8) 



and similarly for k 2 . 

Fig. |3Jc) and (d) shows an endoscope effect wherein 
a point source located close to a tilted array leads to 
focussing and a strong localised beam. The behaviour 
near point M of the Brillouin zone is responsible for the 
focussing effects of Fig. Q]and Fig. EJc). Here the clos- 
est standing wave frequency is at the point M and un- 
equal Ta coefficients that have opposite signs. This leads 
to very strong anisotropy within the effective material 
directed along characteristics; the effective material be- 
comes hyperbolic rather than elliptic. By considering 
only a portion of the Brillouin zone we must bear in mind 
that this asymptotic solution is also valid at N and there 
is a second / equation with the Tn and T%% interchanged 
that is used. The asymptotic PDE for HFH is Eq. © 
with Tn = -9.649 and T 22 = 4.7149. tt 2 is determined 
using the the frequencies stated in Figs [1] and EJc) with 
Eq. ([7|). The opposite sign of coefficients Tn and T22 
yield characteristic type of solutions in order to obtain 
an endoscope effect. 

Last, but not least, the triple crossings in the disper- 
sion diagram Fig. [2] comprise Dirac cones with a flat 
mode passing through the vertex. This is highly inter- 
esting given that Dirac cones are normally limited to 
graphene-like hexagonal structures [17[ ; the current situ- 
ation involves a square lattice is akin to the Dirac cones 
for photonic crystals recently described in [23[ . One 
can use the properties of Dirac cones to reduce the scat- 
tering of a clamped obstacle in a pinned PC, as demon- 
strated in fig. H£a),(b). This behaviour is also captured 
by HFH asymptotics: In this particular case three cou- 
pled equations emerge with variables for i = 1,2,3. 
The system decouples to yield the same governing equa- 
tion for all three functions /i as, 



87T fo^tXi - ^1/0 = 0, 



where Oi is now the first order correction term for the 
frequency The numerics using HFH are also shown 
in Fig. @|c). 

In conclusion, a simplified model of elasticity via the 
thin-plate equation allows for analytical and numerical 
studies opening interesting possibilities in the design of 
flat lens, directive antenna, endoscope, shielding and a 
cloak for flexural waves. Moreover, such designs could 
be scaled up in order to achieve some control of seismic 
surface waves @, [? ], indeed this work is motivated by 
a requirement at CERN to control elastic waves created 



by thermal shock. We also show, for the first time, that 
the computations using the recently developed HFH are 
capable of capturing the fundamental wave propagation 
features of these various possibilities, including counter- 
intuitive physics of Dirac points [23| . 
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